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The FORTE (Fast On-orbit Recording of Transient Events) satellite records bursts of electro- 
magnetic waves arising from near the Earth's surface in the radio frequency (RF) range of 30 to 
300 MHz with a dual polarization antenna. We investigate the possible RF signature of ultra-high 
energy cosmic-ray particles in the form of coherent Cherenkov radiation from cascades in ice. We 
calculate the sensitivity of the FORTE satellite to ultra-high energy neutrino (UHE u) fluxes at 
different energies beyond the Greisen-Zatsepin-Kuzmin (GZK) cutoff. Some constraints on super- 
symmetry model parameters are also estimated due to the limits that FORTE sets on the UHE 
neutralino flux. The FORTE database consists of over 4 million recorded events to date, including 
in principle some events associated with UHE v. We search for candidate FORTE events in the 
period from September 1997 to December 1999. The candidate production mechanism is via co- 
herent VHF radiation from a UHE shower in the Greenland ice sheet. We demonstrate a high 
efficiency for selection against lightning and anthropogenic backgrounds. A single candidate out of 
several thousand raw triggers survives all cuts, and we set limits on the corresponding particle fluxes 
assuming this event represents our background level. 



I. INTRODUCTION 

Detection and modeling of the highest energy cosmic 
rays and the neutrinos which are almost certain to ac- 
company them represents one of the most challenging 
problems of modern physics. To date a couple of tens 
of cosmic ray events, presumably protons, have been ob- 
served with energies in excess of 10^" eV. The origin of 
this flux represents a puzzle since above ~ 5 x 10^^ eV the 
cosmic ray flux is expected to be reduced due to Greisen- 
Zatsepin-Kuzmin (GZK) mechanism. The primary 
particles inevitably generate ultra-high energy neutrinos 
(UHE i^) in cosmic beam dumps. Weakly interacting 
neutrinos, unlike gamma photons or protons, can reach 
us from distant sources and therefore are a possible in- 
valuable instrument of high-energy astrophysics. 

Above ^ 10^^ eV, charged cosmic rays are no longer 
magnetically confined to our galaxy. This implies 
that particles above this energy detected at earth are 
very likely to be produced in extragalactic astrophysical 
sources. Furthermore, the existence of particles up to, 
and possibly beyond the ~ lO"'^^'^ eV endpoint of the al- 
lowed energy spectrum for propagation over cluster-scale 
distances suggests that there is no certain cutoff in the 
source energy spectra. In fact, we are unlikely to learn 
about the endpoint in the source energy spectra via ei- 
ther charged particles or photons: the universe is largely 
opaque to all such traditional messengers. 

There are no clear physics constraints on source energy 
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spectra until one reaches Grand-Unified Theory (GUT) 
scale energies near 10^^ eV where particle physics will 
dominate the production mechanisms. If we assume that 
measurements at ~10^° eV represent mainly extragalac- 
tic sources, then there are virtually no bounds on the 
intervening three or four decades of energy, except what 
can be derived indirectly from upper limits at lower en- 
ergies. For these reasons, measurements which are sensi- 
tive to ultra-high energy neutrinos are of particular im- 
portance in understanding the ultimate limits of both 
particle acceleration in astrophysical zevatrons (acceler- 
ators to ~ ZeV = 10^^ eV and higher energies) and top- 
down decay of exotic forms of cosmic energy. Even at 
energies approaching the GUT scale, the universe is still 
largely transparent to neutrinos. Also, the combination 
of a slowly increasing cross-section combined with the 
large energy deposited per interaction makes such neu- 
trino events much more detectable than at lower energies. 

The chief problem in UHE v detection arises not from 
the character of the events, but from their extreme rar- 
ity. For neutrinos with fluxes comparable to the extrap- 
olated ultra-high energy cosmic ray flux at 10^^ eV, a 
volumetric aperture of order 10^ km'^ sr (water equiva- 
lent) is necessary to begin to achieve useful sensitivity. 
Such large volumes appear to exclude embedded detec- 
tors such as AMANDA or IceCube Q , which are effec- 
tive at much lower neutrino energies. Balloon-based 
or space-based 0, systems appear to be the only viable 
approaches currently being implemented. 

The most promising new detection methods appear 
to be those which exploit the coherent radio Cherenkov 
emission from neutrino-induced electromagnetic cas- 
cades, flrst predicted in the 1960's by Askaryan [H, |||, 
and conflrmed more recently in a series of accelerator ex- 
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periments Above several PeV (lO""^^ eV) of cas- 

cade energy, radio emission from the Askaryan process 
dominates all other forms of secondary emission from a 
shower. At ZeV energies, the coherent radio emission 
produces pulsed electric field strengths that are in prin- 
ciple detectable even from the lunar surface 

These predictions, combined the strong experimental 
support afforded by accelerator measurements, are the 
basis for our efforts to search existing radio-frequency 
data from the Fast On-orbit Recorder of Transient Events 
(FORTE) satellite for candidate ZeV neutrino events. 
Here we report the first results of this search, based 
on analysis of several days of satellite lifetime over the 
Greenland ice sheet. 



II. DETECTOR CHARACTERISTICS 

The experiment described in this paper is based on 
detection of electromagnetic emission generated in the 
Greenland ice sheet by the FORTE satellite. 

The FORTE satelhte JJJ was launched on August 29, 
1997 into a 70° inclination, nearly circular orbit at an 
altitude of 800 km (corresponding to a field of view of 
~ 27° arc distance). The satellite carries two broad- 
band radio-frequency (RF) log-periodic dipole array an- 
tennas (LPA) that are orthogonal to each other and are 
mounted on the same boom pointing in the nadir direc- 
tion. The antennas are connected to two radio receivers 
of 22 MHz bandwidth and center frequency tunable in 
20-300 MHz range. Beside RF receivers, the satellite 
carries an Optical Lightning System (OLS) consisting of 
a charge-coupled device (CCD) imager and a fast broad- 
band photometer. Although for this paper we do not 
report analysis of optical data, in more detailed studies 
the optical instrument data can be used to correlate RF 
and optical emissions [l5l |. 

The satellite recording system is triggered by a subset 
of 8 triggering subbands which are spaced at 2.5 MHz 
separations and are 1 MHz wide. The signal has to rise 
14-20 dB above the noise to trigger. Typically, a trig- 
ger in 5 out of 8 subbands is required. The triggering 
level and algorithm can be programmed from the ground 
station. Multiple channels are needed for triggering be- 
cause of anthropogenic noise, such as TV and FM radio 
stations and radars, which produce emission in narrow 
bands which can coincide with a trigger subband. Af- 
ter the trigger, the RF data is digitized in a 12-bit Data 
Acquisition System (DAS) at 50 Msamples/s, and the 
typical record length is 0.4 ms. The FORTE database 
consists of over 4 million events recorded in the period 
from September 1997 to December 1999. 

The ice with its RF refraction coefficient oin fa 1.8 
and Cherenkov angle of 9c ~ 55.8° and relatively low 
electromagnetic wave losses in the radio frequency range 
is a good medium for exploiting the Askaryan effect for 
shower detection. The biggest contiguous ice volume on 
Earth (the Antarctic) is unfortunately not available to 



FORTE satellite because of its orbit inclination of 70°. 
The next biggest contiguous ice volume is the Green- 
land ice sheet. Its area is 1.8 x 10^ km^, and the depth 
is ~3 km at the peak. However, the available depth is 
limited by RF losses jT^ to 1 km. Thus the volume of 
Greenland ice observed from orbit is w 1.8 x 10® km^^. 



III. CHERENKOV RADIO EMISSION FROM 
PARTICLE SHOWERS 

We define the electric field pulse spectrum as E(ct;) = 
2 E(i)e*"*(it. An empirical formula for E(aj) from 
an electromagnetic shower in ice was obtained by Zas et 
al. d: 
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where R is the distance to the observation point in ice, 
i^^showcr is the shower energy, assumed to be <1 PeV, / 
is the electromagnetic wave frequency, /o = 500 MHz, 
Oc = 55.8° is the Cherenkov angle, and AO = 2.4°/o//. 

The FORTE detector triggers whenever the amplitude 
of the electric field after a narrow-band (1 MHz) fil- 
ter exceeds a set threshold (in several channels). Since 
|E(a;)| varies slowly enough in this band, the peak value 
of E on the filter output equals |E(w)| A/, where A/ and 
/ — ti)/(27r) are correspondingly the bandwidth and the 
central frequency of the filter. 

For low frequencies the emission cone is broad 
{A0 ~ radian) and the empirical formula ^ is not very 
accurate. Instead, we make an analytical estimate for 
emission in Appendix IXI and get 
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where k = i'Knf /c and n — 1.8 is the ice refrac- 
tion coefficient, and ^ = 1. At fcL ^ 1 this formula 
matches the empirical formula for L 1.5 m and 
Q w (5showor/l TeV) X 5.5 X 10^^'' C, which agrees with 
the results of shower simulations in Figures 1 , 2 in . 
Equation ^ also matches the numerical result for the ra- 
diation pattern at 10 MHz (Figures 11, 12 in [ll|) better 
than the empirical formula ^ which becomes inaccurate 
at low frequencies. 

Note that we cannot use equations Q for elec- 
tromagnetic showers started by particles of high ener- 
gies because of significant elongation due to Landau- 
Pomeranchuk-Migdal (LPM) effect [H H HI ll^ . Ac- 
cording to [23, the LPM effect is important for particle 
energies S > ^lpm, where Slpm = 2.4 PeV. According 
to Appendix^ for electromagnetic showers with start- 
ing energy Sq ^ ^lpm, the electric field at the exact 
Cherenkov angle is still approximately given by (Q, while 
the width of Cherenkov cone is reduced to 

A0LPM~O.9° ^ I ^ (3) 
y/£o/l EeV / 
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where /o = 500 MHz. For this reason, pure electromag- 
netic showers (primarily from v^, or charged current in- 
teractions) make a negligible contribution to the FORTE 
sensitivity because, as we will see, the energy threshold 
is >1 EeV. 

After the interaction with a nucleon, the UHE v energy 
goes into leptonic and hadronic parts. In the case of an 
electron (anti-)neutrino i^e iye)-, the electron created in 
the charged current interaction starts an electromagnetic 
shower, which is however very elongated due to LPM 
effect making it virtually undetectable. In the case of a 
v^i (f'/j) or Vt (vr) the lepton does not start a shower. 
The muons and tau leptons deposit their energy due to 
electromagnetic Q and photonuclear interactions. 
However, the portions in which the energy is deposited 
are not enough to be observed by FORTE. Also, the tau 
lepton may decay at a long distance from its creation 
point after its energy is reduced to ^1 EeV (10^* eV). 
Thus, tau lepton decay is also unobservable by FORTE. 

On the basis of these arguments, we only consider the 
neutrino-initiated hadronic shower, and we do not dis- 
tinguish between neutrino flavors. Most of the hadronic 
energy converts in the end into electromagnetic, thus 
Cherenkov emission as a result of Askaryan effect can 
be used for its detection. Moreover, the LPM effect does 
not produce significant shower elongation, as shown by 
Alvarez-Muhiz and Zas j26| , since most of the ttq particles 
which decay into photons instead of interacting have their 
energy reduced below the LPM level. Thus, the value of 
L « 1.5 m used in Q does not change appreciably. This 
is in accordance with results of Monte Carlo calculations 
of '26^1 , showing that the Cherenkov cone narrows by only 
~30% when hadronic shower energy ranges from 1 TeV 
to 10 EeV. 

The hadronic shower energy which should be used in 
or is £showcr = v£-v, where Ey is the neutrino 
energy, and y is the fraction of energy going into the 
hadronic shower. The theoretical value for UHE v is 
(y)«0.2 H. 

Because the long wavelength observed by FORTE 
(~6 m in ice) are far greater than the Moliere (transverse) 
radius of the showers involved (^11.5 cm), the trans- 
verse structure of the shower (expressed by the factor 
(1-1-0. 4(///o)^)"^ in (UJ) is neglected in our analysis. The 
same argument applies to transverse non-uniformities of 
hadronic showers. 



IV. SEARCH FOR RELEVANT SIGNALS IN 
THE FORTE DATABASE 

A. Geographic location of FORTE events 

The geographic location of the signal source can be 
determined using the dispersion of the short electromag- 
netic pulse in HF range going through ionosphere. Two 
important parameters used in geolocation of the source 
can be determined from the data from a single FORTE 



antenna. The first parameter is the total electron con- 
tent (TEC) along the line-of-site between the source and 
the satellite. It is proportional to the group time delay, 
which has frequency dependence |28l |. The second 
parameter is determined from the Faraday rotation of a 
linearly polarized signal 29] , due to birefringence in mag- 
netoactive ionospheric plasma. The Faraday rotation fre- 
quency turns out to be equal to the "parallel" electron 
gyrofrequency J\\^ce = eB\\/me = fee cos 9, where is the 
angle between the geomagnetic field B and the ray trajec- 
tory at the intersection with ionosphere. Both frequency- 
dependent delay and frequency splitting due to Faraday 
rotation are well seen in Figure ^ The Cherenkov radio 
emission is expected to be completely band-limited and 
linearly polarized, which enables us to make use of the 
second parameter for geographic location. 

We calculate the probability distribution of the source 
location using Bayesian formula: 

p({A,0}|TEC,/||,,e) = Cp(TEC|{A,0})p(/||.,,|{A,0}) 

(4) 

where {A, (j)} are the latitude and longitude of the source, 
p(TEC|{A, (/)}) and p(/|| ce|{A, 0}) are conditional proba- 
bility distributions for the measured parameters given the 
location of the source, and C is a normalization constant. 
Here we assume that the measurements of parameters are 
independent, and the a priori distribution of the source 
location is uniform in the field of view of the satellite. 

To estimate TEC between the source and the satellite 
(for given locations of both), we use the Chiu ionosphere 
model 30], adapted to IDL from the FORTRAN source 
code found at NASA ionospheric models web site. This 
model gives electron density as a function of altitude for 
given geographic and geomagnetic coordinates, time of 
year, time of day and sunspot number. By integrating 
it over altitudes, we find the vertical TEC. To convert it 
to TEC along the line-of-sight, we must divide it by the 
cosine of the angle with the vertical. Due to the curva- 
ture of the Earth, this angle is not constant along the 
line-of-sight, and we make an approximation of taking 
this angle at the point where the line-of-sight intersects 
the maximum of the ionosphere (-F-layer) , at altitude of 
^300 km. The Chiu model, due to simplifying assump- 
tions, does not account for stochastic day-to-day variabil- 
ity of the vertical electron content. The standard devia- 
tion can be as large as 20-25% from the monthly average 
conditions (3 It p. 10-91]. Thus, we assume a Gaussian 
probability distribution for p(TEC|{A, (/>}) with the cen- 
ter value calculated using Chiu model and variance of 
25%. 

The geomagnetic field is estimated from a simple dipole 
model py, pp. 4-3, 4-25]. The error is assumed to be 10% 
according to the estimates for experimental determina- 
tion of /ii^ce from the RF waveform, in Figure 6 of [2^ . 
However, this uncertainty can be greater for signals that 
are only partially linearly polarized. Again, we use a 
Gaussian distribution forp(/|| ,^g]{A, </>}) with correspond- 
ing central value and the standard deviation of 10%. 
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FIG. 1: An example of a highly-polarized impulsive event detected by FORTE. The confidence levels for event location 67% 
and 90% determined using equation Q S'l^s shown in the map. The spectrogram shows the dispersion of the pulse in the 
ionosphere and splitting due to Faraday rotation in geomagnetic field. The horizontal lines are due to anthropogenic noise (TV 
and FM radio stations). 



B. Background rejection 



FORTE SATELLITE SENSITIVITY TO 
NEUTRINO FLUX 



The pulse generated by a UHE v shower in ice is 
expected to be highly polarized and essentially band- 
limited up to a few GHz. In these aspects, it is simi- 
lar to the electromagnetic emission from the "steps" in a 
stepped-leader lightning "s^]. However, the pulses corre- 
sponding to lightning steps are accompanied by similar 
neighbors before and after, within a time interval from a 
fraction of a ms to ^0.5 s. The signal grouping can thus 
be used to distinguish UHE ly signatures from most light- 
ning events. Also, the lightning activity must be present, 
which is extremely rare in areas of the Earth covered by 
ice, and thus can be excluded using the method described 
in Section HVH 

There is a special type of intracloud lightning which 
produces isolated events which are called compact intr- 
acloud discharges (CID) (33|- However, these events are 
usually randomly polarized and have several-/zs pulse du- 
rations |3^. 

Another rejection method uses the fact that the light- 
ning discharges occur above ground, and therefore there 
is a large probability for FORTE to detect also the signal 
reflected from the ground. This phenomenon is known as 
Trans-Ionospheric Pulse Pairs (TIPPs) [Ulllllll. The 
presence of a second pulse, therefore, excludes the possi- 
bility of the signal to be a UHE signature. An example 
of a TIPP event in FORTE data is shown in Figure |21 



In this section, we estimate the upper limit of the dif- 
ferential flux of UHE ly (all flavors) depending on the 
number of triggers on the relevant events in the satellite 
lifetime. The limits are set on the sum of all flavors since 
this experiment cannot distinguish between them. 

The typical natural background noise level in FORTE 
data is -lO-^^^iO-" (V/m)2/MHz (as can be seen, e.g., 
from spectrograms in this paper's Figures). In a typ- 
ical 1-MHz trigger subband this corresponds to RMS 
value of 1-3 fiV/m. The trigger level is set 14-20 dB 
above the noise, giving the ability to trigger on impul- 
sive signals with frequency domain values of 5-30 /LtV/m 
in each 1 MHz trigger subband. For flux limits cal- 
culations in this section, we use the threshold value of 
30 ^Vm^i MHz^i at / = 38 MHz, the central frequency 
of the low FORTE band. 



A. Relation between FORTE sensitivity and the 
limits on UHE i/ flux 

Let X{£) be the theoretical number of triggers of 
FORTE satellite in its full lifetime assuming a unit mo- 
noenergetic neutrino flux for different energies £. We will 
call this function the sensitivity of the FORTE satellite 
to UHE ly flux. Then the expected number of triggers is 
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FIG. 2: A TIPP (pulse pair) event with probable origin location in Greenland but rejected in our analysis. Notations are the 
same as in Figure 



where ^{E) is the differential neutrino flux (per unit area, 
time and sohd angle). The number of detected events is 
a Poisson distributed random number with expectation 
s. If no events are detected (null result) we can set a 
limit on s: 

s < Sup = - log a 

where 1 — a is the confidence level. For 90% confidence 
level, for example, we have 

J X{£)<^>{£)d£ < Sup w 2.3 (5) 

In general, for n events, the limiting value of s is deter- 
mined from 

a = (3(n + 1, Sup) (6) 
where Q(n -\- 1, Sup) = -4^ is the regularized in- 

I x'^e ^ dx 
Jo 

complete gamma function. For example, at a confidence 
level of 90% and n = 1 we get Sup ~ 3.89. 

From the limit set on integrated flux jSj) one can try 
to construct a limit on a differential flux. However, such 
a limit can be evaded for differential flux models that 
are anomalous, for example, if there are very narrow 
emission lines in the spectrum. Nevertheless, for reason- 
able assumptions such as a smooth and continuous model 
spectrum, the implied limits are model- independent. As 
shown in Appendix |C1 if we assume that the spectrum is 
sufficiently smooth, i.e. does not have any sharp peaks 
(the peaks have widths at least of the order of the central 
energy of a peak) , then we can assert that 



Since this limit does not assume any particular model 
(except that it is sufficiently smooth) , we will define it as 
the model-independent limit. 



B. Sensitivity calculation 

The sensitivity \(£^), which was introduced in the pre- 
vious subsection, is calculated using the specific aperture 
R, which we define as the trigger rate for unit monoener- 
getic neutrino flux when neutrinos interact under a unit 
area of ice. The specific aperture R is integrated over 
the visible ice area averaged over time (many satellite 
passes) and multiplied by the entire satellite lifetime in 
orbit T (while the two radio receivers of 22 MHz band 
were in order) and the fraction of time in trigger mode 
(duty cycle) D: 

\'^'4vi5ibic / time 

The hfetime of FORTE T is from September 1997 tiU 
December 1999. The time when the emission from the 
Greenland ice sheet is visible by the satellite can be es- 
timated by calculating the total time spent by FORTE 
inside a circle of radius 20° arc distance from the approx- 
imate geographic center of Greenland, at 70°N, 40°W. 
This time is estimated to be ~38 days. The duty cy- 
cle D is estimated to be 6% by calculating the run time 
fraction within this circle, making the effective time of 
observation to be ~3 days. 

The specific aperture i? is a function of the position of 
the neutrino-generated shower which is given by the arc 
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FIG. 3: Configuration for sensitivity calculations. 

distance to the shower location 9^ and its azimuthal an- 
gle (f)a (see Figure E)). Note that the visible area Avisibio 
depends on the satellite location. Averaging over time 
in our calculations is replaced by averaging over satellite 
position over the Earth surface with a weight correspond- 
ing to the fraction of time the satellite spends in a given 
point of the globe. 

To calculate the specific aperture R, we must calcu- 
late how many neutrinos interact in ice, depending on 
energy dip angle adip (the angle of v velocity below 
the horizon at the interaction point) and the depth of 
interaction z. We use the theoretical neutrino- nuclcon 
interaction cross-section (T,y AT (£j/) The neutrino flux 
from below is greatly reduced by interactions in the Earth 
volume, and most detectable interactions are from hor- 
izontal and down-going events. Here we have to note 
that this statement does not apply to neutralino inter- 
actions since neutralino-nucleon interaction cross-section 
a^M can be much smaller than neutrino cross-section at 
a similar energy js^. 

The number of interactions in ice is characterized by 



function F{Vl. z\£^) defined in the following way: 

^F{n,z;£^)dn dV = Nnuc(7^Ni£u)d<P dV 

Here fl = (adip, 0s) is the neutrino velocity direction 
(see Figure 01), z is the interaction depth, N^uc is the 
nucleon number density and d^ = ^e~'^d^} is the flux 
in solid angle element dfl = sinadipdadip'i'/'s attenuated 
by neutrino absorption in a layer of optical thickness 

T(adip,z) = / N^^cCfyNdl. 

The specific aperture R is found using the fraction of 
interacting neutrinos which produce field i?ant > ^'th at 
the satellite antennas: 

RiOsAa^E.) = J F{n,z;£,) 

^ e[E,nt{y£., ^, 9s, (j)a) - EthMy) dy^ dfl dz 

where Q is the step function, and p{y) = {1/ a){dij / dy) 
is the probability distribution function for the kine- 
matic parameter y = £h_ad/£u, the fraction of energy 
going into the hadronic shower. The p{y) dependence 
is taken from for the highest energy considered in 
that paper, E^, = 10^^ GeV. Since Es,ntiy£u,^,9s,(t)a) = 
yEaati^f, ^,9s,4>a), we can integrate over y immediately 
and get 

R{9s,ct>a-£.) = j F{n,z-£,) 

.'-^^(i?ant(jj^,^.>.) 

where Fy{y) — p(y')dy' is the cumulative distribution 
function of y. 

Given a neutrino interacting at dip angle adip and 
depth z in ice, we can calculate the field detected at the 
satellite. First, the emitted field is given by Cherenkov 
emission formula |(2J). Then the field is attenuated in ice, 
refracted at the ice surface and detected at the satellite 
taking into account the directional antenna response. In 
the subsections of Appendix^ we give the details of cal- 
culations for these steps. In Figure 01 we show the peak 
electric field at the satellite altitude (800 km) at the cen- 
tral frequency of the low FORTE band (/ = 38 MHz). 
As we see, even for a dip angle of 20° there is a signifi- 
cant emission upward from the shower due to the width 
of Cherenkov cone. Although the field plotted in Fig- 
ure 01 does not include the antenna response, the plots 
can give some idea about where the satellite has to be to 
see the signal at a given threshold, the typical value of 
the threshold being 30 /j,Vm~^ MHz~^. 

VI. RESULTS 
A. Event search results 

We searched for events recorded while FORTE was in- 
side a circle of radius of 20° with a center at 70°N, 40° W, 
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(a) PeakE, ^Vm^MHz^forE^g^ = 2x10^^eV (b) Peak E, m'^MHz^for E^3^ = 2x10^^eV 




FIG. 4: The peak electric field at the satellite altitude -Esat, emitted in a Cherenkov process by a hadronic shower of energy 
2 X 10^^ eV. These plots do not include the antenna response. The shower is in the center of the circle, and is directed to the 
right. Emission is shown for interaction at depth z = 100 m and dip angles of (a) Odip = 0° and (b) Odip = 20°. The dashed 
circles represent the arc distance in degrees (1° « 111 km on the Earth surface). 



in time period from the start of FORTE in September 
1997 to December 1999, when both 22-MHz-bandwidth 
receivers were lost 39j. We estimate that the satellite 
spent a total of 38 days inside this circle, with at least 
~6% of it being the time in trigger mode. We found a 
total of 2523 events. From these, only 77 are highly po- 
larized. These 77 events can be geolocated using both pa- 
rameters described in subsection II VAI Of these, only 16 
events have intersection of the 90% confidence level with 
Greenland's ice sheet. Out of the remaining 16 events, 11 
are rejected for being TIPP events, i.e. pulse pairs with 
ground reflections that indicate that the origin locations 
are above ground. An example of a rejected TIPP event 
is shown in Figure |21 Two more events were rejected 
because of the presence of a precursor before the pulse, 
which is characteristic of certain type of lightning 
and cannot be present in a neutrino shower signal. An 
example of such event was shown in Figure ^ 

Out of the remaining three events, one (shown in Fig- 
ure is rejected for its long duration (>10 /is), since 
the Cherenkov pulse is expected to be only ~1 ns long 
(the time resolution of the FORTE detector limits this 
to -^20 ns) . The remaining two events are shown in Fig- 
ures El and [7| The first of them has close neighbor events 
(at —1.4 ms and -1-0.7 ms), which makes it a probable part 
of a stepped-leader process in a lightning. The neighbors 



of the second event are not very close (at —0.27 s and 
-1-5.55 s), but it still can be a lightning event. The recent 
analysis by the FORTE team has shown that the 
lightning events are likely to have neighbors in ±0.5 s in- 
terval, with the most probable separation of ±0.01 s, and 
the accidental coincidence rate is ~0.9 per second. This 
value of the accidental coincidence rate makes the candi- 
date event shown in Figure |7| indistinguishable from an 
isolated lightning discharge. The analysis of this event 
continues. 



B. Flux limits set by FORTE 

In TableUwe give the calculated FORTE neutrino flux 
sensitivity values XiEu). On the basis of these values one 
can set a limit on any model flux ^(S-^) using numerical 
integration in equation jSJl, with s^p — 3.89 since we 
have one uncertain event as our background noise. In 
FigurelHlwe plot the model-independent flux limits set by 
equation |(7J) on the basis of these data (again with Sup — 
3.89). In the same Figure we also show the comparison 
of calculated flux limits with predicted neutrino fluxes 
from various sources. Some of the sources of super-GZK 
neutrinos are reviewed, e.g. in [4]| . As we see, the flux 
limits set by FORTE observations of the Greenland ice 
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FIG. 5: The confidence levels of 67% and 90% for geographic location and a spectrogram of an example event. This event 
cannot be generated by a neutrino due to its long duration (^10 /is). 
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FIG. 6: The confidence levels of 67% and 90% for geographic location and a spectrogram of an example event. A short 
horizontal streak in the spectrogram is due to anthropogenic noise (a radar). This event needs further consideration for being 
rejected as neutrino-generated. However, its nearest neighbor events were found at —1.4 ms and +0.7 ms, which makes it a 
probable lightning event. 



sheet can reject some regions of parameters of the Z burst 
model milllillillMlil. 

Note that differential fluxes in some models can be even 
smoother and wider in energies than assumed for deriva- 



tion of ^ . Figure |2| shows what would be the hmits on 
a class of models with power-law energy dependence of 
the flux, oc £^". A similar analysis using power- 

law models was performed in the past, e.g., for RICE 
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FIG. 7: The confidence levels of 67% and 90% for geographic location and a spectrogram of an example event. Its nearest 
neighbors were found at —0.27 s and +5.55 s. This event needs further consideration for being rejected as neutrino-generated. 



TABLE I: The FORTE sensitivity X{£) (cm^ ssr) to neutrino 
flux (any neutrino flavor) and neutralino flux (for different 
neutralino-nucleon cross-sections) . 
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detector [ssf. 

As another result of present research, we set limits on 
the flux of neutralinos, weakly interacting particles pre- 
dicted by the Minimal Supersymmetric Standard Model 
(MSSM). The calculations are performed in the same 
way as for setting the neutrino flux limits with a few dif- 
ferences. First, the cross-section for nucleon interaction 
a^N is different, and is expected to be in the range from 
~ (l/100)(TyAr to ^ (7^N [33. Second, all of the energy 
goes into the hadronic shower. The results are presented 
also in Table H In Figure ITUl we compare the predicted 
neutralino fluxes 38] with the model-independent flux 
limits set by FORTE. We see that for neutralino-nucleon 



cross-sections in the range a^N O.la^N, strong lim- 
its are set on the model predicting neutralino flux from 
decay of heavy X particles with Mx = 2 x 10^^ eV, espe- 
cially if the sources are homogeneously distributed. Note 
that for a case of a different mass, Mx — 2 x 10^^ eV, 
which is also considered in , we are unable to set any 
limits since all decay products arc below FORTE thresh- 
old. 

In Tabled we apply the more robust model-dependent 
limit to these models to get the confidence level of rejec- 
tion according to equation ((^Jl, for n — I (one uncertain 
event). We vary cross-section a^x, take different decay 
scenarios and X distributions. As one can see from this 
Table, several models are rejected with very high confi- 
dence. 



VII. DISCUSSION 

The limits shown in Figures |H1 and ^1 represent to 
our knowledge the first direct experimental limits on the 
fluxes of neutrinos and other weakly interacting particles 
in this energy range. The fact that the first such limits 
already have constrained several proposed models is an 
indication of the power of the radio detection techniques, 
but the scarcity of other limits in this regime also sug- 
gests that they be accepted with caution. Here we discuss 
briefly some of the potential issues with these constraints. 

At the energies to which FORTE is sensitive, the 
energy of the pulsed coherent radio emission can be- 
come one of the dominant energy-loss mechanisms for 
the shower. This implies that the radiation reaction of 



10 



TABLE II: The confidence levels for rejecting the models of neutralino production by heavy X particle decay with AIx = 
2 X lO'^^ eV fs^. We only show models with the rejection confidence level >50%. The scenarios of X decay are the same as in 
Figure nUl The variable a is defined as 1 — CL. 
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FIG. 8: The estimated limits on UHE v flux detectable by 
FORTE using Greenland ice sheet. The limit is compared to 
predicted neutrino fluxes from various sources: GRB: gamma 
ray bursts ^481; AGN: active galactic nuclei .491 : GZK: neu- 
trinos produced in GZK mechanism from cosmic rays fsof : 
TD: topological defect model (non-SUSY) |5l|: Z burst mod- 
els are: the slanted box bv |44l. the diamonds showing the 
la level errors are by |4^ l52(l . and the arrow is by [isl de- 
fault value of parameters]. The limits from other experiments 
are also shown: AMANDA [13, RICE [H (determined us- 
ing (|7|| from the effective volume and time of observation), 
GLUE [T3. 17^ . The limits from all experiments show the 
limits on the combined flux of neutrinos of all fiavors (as- 
sumed mixed in equal amounts), except AMANDA which is 
only sensitive to i/f^ (^p.)- 



the shower to the pulse could lead to modification of the 
shower development, and consequently some change in 
the radiation parameters. We have not attempted to cor- 
rect for this effect in our analysis, but we note that it is 
probably not important below ~ lO^"* eV T^. Above this 
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FIG. 9: Upper curve: the differential fiux limit of equation 
||7|l; lower set of lines: the limits set by equation lICH assum- 
ing power law flux shapes K[£; a) = £~" for a in the interval 
from 1.6 to 6.0 with a step of 0.4. 



energy we expect that the shower radiation pattern might 
spread to some degree, depending on the foreshortening 
of the cascade length due to radiation reaction decelera- 
tion. 

We have noted that our limits extend into the mass 
scale for GUT particles. However, it is important to note 
that the center-of-momentum energies for interactions of 
these neutrinos on any other standard particles are of 
order 10 PeV or less. Our results thus depend on a <3 
order of magnitude extrapolation of the standard model 
neutrino cross sections from the current highest energy 
estimate from accelerators at ~30 TeV [Sa, |53, llll , over 
an energy regime where the cross sections grow only log- 
arithmically with energy. For this reason we do not ex- 
pect that the energy scale itself is good cause to doubt 
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FIG. 10: The limits on neutralino fluxes set by FORTE ob- 
servations of the Greenland ice sheet, for one detected event 
and different assumptions about neutralino-nucleon cross- 
sections. Shown also are predicted neutralino fluxes l3a l 
for decay of superheavy particles of Mx = 2 x 10 eV. 
The lower four curves are for X decays in the halo of the 
Galaxy, with the primary decay into: (1) quark-|-antiquark 
(solid); (2) quark-|-squark (dot-dash); (3) 5(7(2) doublet 
lepton-|-slepton (dots); (4) 5 quark-|-5 squark (short dashes). 
The upper curve (long dashes) is for homogeneous X distri- 
bution (in which case the flux is enhanced by a factor of 15 
compared to a "galactic" distribution), decay scenario 3. See 
also Table ini 



scenarios consistent with FORTE data. 

The strong FORTE constraints on the neutralino pro- 
duction model [111 from heavy X particles, sets a joint 
limit on (1) X particle distribution and mass (2) X par- 
ticle decay channels and (3) neutralino-nucleon cross- 
section. Since some of X particle decay scenarios (at 
the mass 2 x 10^^ eV) are strongly rejected, as shown 
in Figure ^| and Table ^ this can give an insight into 
the possible nature of such particles and the physics at 
supersymmetric grand-unification scale. Although refer- 
ence |38j does not consider models with masses interme- 
diate between 2 x 10^^ and 2 x 10^^ eV, we expect that 
forte's sensitivity, which extends an order of magni- 
tude or more below the region where the model at 2 x 10^^ 
is constrained, will also limit neutralinos of masses down 
to ^ 10^^ eV, particularly if the cross sections approach 

VIII. CONCLUSIONS 

We have performed a search for radio frequency sig- 
natures of ultra-high energy neutrinos originating from 
coherent Cherenkov emission from cascades in the Green- 
land ice sheet, observed with the FORTE satellite over 
an ~2 year period. In ~3 days of net exposure, a sin- 
gle candidate, presumed to be background, survives the 
analysis, and we set the first experimental limits on neu- 
trino fluxes in the W^-W^ eV energy region. These 
limits constrain the available parameter space for the Z 
burst model. In addition we constrain several variations 
of a model which involves light super-summetric particles 
(neutralinos) at these energies, particularly those with in- 
teraction cross-sections approaching those of neutrinos. 



the values for the flux limits. 

At 10 PeV center-of-momentum energies, interactions 
of the primary neutrinos will exceed in CM energy those 
of any observed ultra-high energy cosmic rays, and could 
therefore lead to production of new heavy particles in 
the showers themselves. Such interactions could include 
channels in which most of the energy goes into an un- 
observable particle, or a particle with interactions much 
weaker than neutrinos. In the absence of any specific 
proposals for models and interactions, we can only note 
that such behavior can evade our limit, but could lead to 
other observable secondary particles with different angu- 
lar distributions that we have not considered. 

Since FORTE threshold for detection of a weakly in- 
teracting particles (~ 10^^ eV) is higher than the GZK 
limit (~5 X 10^^ eV), it sets the limits on neutrinos pro- 
ducing super-GZK cosmic rays through resonant interac- 
tion with background neutrinos within ~50 Mpc distance 
from us, the Z burst mechanism [i l li HiHSH i R lir f . 
Although FORTE sets limits on parameters of Z burst 
models, the uncertainties in the models and the mea- 
sured super-GZK cosmic ray flux still make most Z burst 



APPENDIX A: CHERENKOV EMISSION FROM 
ELECTROMAGNETIC SHOWERS 

Let us model the shower as a point charge moving with 
the speed of light. Then the current is given by Jz (r, t) = 
cq{z)b(v — czt) and — Jy ^ 0. The Fourier transform 
is 

Mr,cu) = 2 J J^e'^'dt = 2q{z)S{x)S{y)e'^'^' 

(The factor of 2 is for consistency with our definition of 
E{ll})) The frequency-domain vector potential A satisfies 
Helmholtz equation W'^Az + k^A = — /i/io and A^ — 
Ay — 0, where k = nuj/c and n — -y^e/i- Its solution at 
the observation point R is 

A.(R)=MAio J ^Mr)d''r 

where R' = |R — r|. In the Fraunhofer zone, the standard 
approximation is 

£ ^ f g-»(k-r) 

47ri?' ^ AttR 
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Thus, 

ikR r+oo 

where 6 is the emission angle. 

The magnetic induction is B = V x A. In the far zone 
B = |B| = kAz sin0 and the electric field \s E ~ cB/n — 
LoAz sinO. 

Let us consider a Gaussian shower profile q[z) — 

Qe ^ where Q is the maximum attained charge ex- 
cess and L is the characteristic shower length. Then 



i£|E(a;)| = ^ — smf 



APPENDIX B: 



-{kLy {cose-l/ny/2 



CHERENKOV EMISSION FROM 
LPM SHOWERS 



According to |23| , LPM effect is important for particle 
energies £ > ^lpm, where ^lpm = 61.5Ao TeV, where 
Xq is the radiation length in cm. The radiation length in 
mass units is 36.1 g/cm^ in water, giving Xq — 39.1 cm 
in ice since picc — 0.924 gcm"'^. Thus, Slpm — 2.4 PeV. 
The increased radiation length for bremsstrahlung and 
the 4/3 of the mean free path for pair production, accord- 
ing to the same paper, are given approximately (within 
-^20%) by Alpm = \/£/£lpmXo, when £ > £lpm- Let 
us model the UHE electromagnetic shower as the initial 
particle gradually losing its energy, which goes into pro- 
duction of usual "small" NKG showers each having an 
initial energy of £lpm- Using this information, we can 
write the energy loss equation: 



LPM 



&t £ > £] 



LPM 



LPM 



where t — z/Xq is the thickness in radiation lengths. 
Solving this equation, we find the number of showers per 
unit length with starting energy f lpm : 



1 



d£ 



.-LPM 



Xodt 



for ti <t <ti + 2tLPM 



where ^lpm = -'^lpm.o/^'^o = V^o/^lpm, and U is 
the depth of the first interaction which can be taken 
U = ^lpm- The number of particles in each "small" sub- 
shower can be described approximately as 



Nsmiz) 



where L — 1.5 m as established in Section IlIII and Zmax 
is the location of the maximum. The maximum number 
of particles is given by p. 23] 



N„ 



37V, 



1 



£ 



^\og{£/£,) £c 



where £c = (610 MeV)/{Z + 1.24) is the critical energy 
[eot- which for water {{Z) = 7.22) is equal to 72.1 MeV. 



At £ = £lpm , we have N„ 



,lpm 



8x10''. We assumed 



that there are equal numbers of electrons, positrons and 
photons. The total number of particle in the shower is 
given by a convolution: 

N{z)= J Psra{z')N,r,^{z~ Z')dz' 

See comparison of this approximate theory and the 
results of Monte Carlo calculations using program 
LPMSHOWER [63 in Figure CH 

The charge excess is estimated to be 



0.2 



Thus q{z) K, 0.2e(2/3)A^(z) where e is the electron 
charge. By the way, &i £ = 1 TcV, we have gmax — 9.6 x 
10"" C, almost the resuh of Section|lIllof 5.5 x lO^i^ C. 

Let us take its Fourier transform, N{p) — 
J N{z)e~'^P^dz. We use the fact that the Fourier trans- 
form of a convolution is just the product of Fourier trans- 
forms. The needed Fourier transforms are 



Psm[P) — e ^LPM 

where a = pXi^pufl, and 

Nsm{z) = e-'P^-'-^TVinax.LPM 



sma 



sma 

a 



cos a 



For our purposes, the exact absolute phase is not impor- 
tant. The electric field is 

i?|E(w)| = ^^og(fc[cos6i- l/?i])/sin6' 

At p = (i.e., Cherenkov angle), we get the maximum 
value 



i?|E(c^)|, 



1.4 X IQ-^- 



£o 



1 TeV /o 



VMHz-i 



which is approximately the same as before. However, the 
width of Cherenkov angle at high energies is determined 
by the extended length of the shower. We can define it 
as the angle at which |E| is reduced by a factor of e~^^^, 
which occurs at a = pAlpm.o ~ 2. We approximate 
p = k{cos{9c + AO) -1/n) fa (27r/n/c) sin^c and get 
the cone width due to LPM effect to be 



where /o = 500 MHz. 



0.9° 



1 



/o 



y/£o/l EeV / 



APPENDIX C: MODEL-INDEPENDENT LIMIT 
ON DIFFERENTIAL FLUX 

As we mentioned in Section from a single equa- 
tion (jSJ one cannot set in general a model-independent 
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FIG. 11: Results of Monte Carlo simulation compared to simple analytical model presented here. 



limit on a differential flux ^{£)- However, after cer- 
tain assumption of smoothness of function it can 
be done |62l |. Let us first consider a model- dependent 
limit on differential flux from a single condition © 
assuming that <&(£) has a certain functional form. Usu- 
ally, it is assumed that = ^oK{£;P) where K is 
a functional shape determined by a set of parameters P. 
Then from ^ it follows that 



APPENDIX D: 



'up 



< max 



J Xi£)K{£;P)d£ 
s^pK{£; P) 



p J X{£')K{£';P)d£' 



(CI) 



It turns out that this equation is valid even when $(f ) 
is a linear combination of functions K{£; P): 



MP)K{£;P)dP 



We can prove it assuming the opposite. If 
$(£) / \{£')K{£'; P)d£' > s^pK{£; P) for aU P, then 

^{£')\{£')d£' = j <^o{P) (^j X{£')K{£';P)d£'^ dP 



J^o{P)K{£;P)dP 



which contradicts our initial assumption (|3J). 

Although in this paper we will not use any concrete 
functions K{£;P), we get a simple formula from 
(|Cip by assuming that K{£;P) = K{£;£o) is a curve 
of width of ^ £o centered at £q and normalized so that 
J K{£;£o)d£ = 1. Then maxg^^ K {£ ; £o) is achieved at 
£o £ and is ^ l/£, and if X{£) is smooth enough, 
/ X{£')K{£';£o)d£ w X{£o), and is w X{£) when the ex- 
pression on the right-hand side of ljCl|> is maximized. So 
we estimate 



< 



~up 



£Xi£) 



Thus, we have shown that a certain region of differen- 
tial fluxes can be rejected on the assumption that they 
are sufficiently smooth functions of energy £. 



ELECTRIC FIELD AT THE 
SATELLITE 



1. Transmission through ice 

The Greenland ice sheet at depths <1000 m has tem- 
peratures from -25° C to -20° C [13, pp. 23-24]. The 
attenuation in ice at frequency / = 35 MHz at these 
temperatures is given by [161 and is «1 dB/100 m. 



2. Refraction 

First, let us find the refraction angle r. Consider satel- 
lite at altitude hgat , and the particle shower occurring at 
arc distance 9s from satellite position. Since the depth 
of the shower is small compared to the satellite altitude, 
we can assume that refraction also occurs at arc distance 

The distance from the satellite to the refraction point 
is found using cosine theorem: 



Rs 



Rl 



{Ra 



2i?e(i?e -I- hsa,t) cos 6s 



where i?^ is the Earth radius. The nadir angle b is found 
from sinfe = R^sinOs/Rsnt- The refraction angle r can 
be found from i?0sinr = (i?^ -I- /isat)sin6, and the in- 
cidence angle from Snell's law, sini = sinr/n. Since 
i,r £ [0, 7r/2], cosi, r — \/ 1 — sin^ i, r. 

After the refraction at the ice surface, E changes ac- 
cording to the Fresnel formulas [63l pp. 281-282] 

E'j_ 2n cos i 

E± n cos i + cos r 
E'^l 2n cos i 



'II 



n cos r -|- cos i 



where i and r are the angles of incidence (from below) and 
refraction, correspondingly, related by the Snell's law, 
nsinz = cosr; E and E' are the incident (below the sur- 
face) and refracted (above the surface) electric field com- 
ponents, E± and E|| are the components perpendicular 
and parallel to the plane of incidence. 
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However, we also need to know how the waves diverge 
to be able to use the expression for i?|E|. Let R be the 
distance from the source to the point at which refrac- 
tion occurs. If we look from above the surface, the waves 
diverge in such a way that they look like coming from 
distance R' below the surface. Then at the satellite, the 
field is determined by relation E^atRsat = E'R', where 
-Rsat ^ R,R' is the distance to the satellite. The in- 
equality is well justified since the shower occurs at depth 
~1 km in ice, while the satellite altitude is 800 km. To 
find R', consider an area element dA of the surface. Then 



dA 



R^dn R'^dfl' 



cost 



COS r 



where dJ7 is the solid angle element at which dA is seen 
from the source point and dfl' gives divergence of rays 
emanating from dA above the surface. These solid angle 
elements are dfl — sin i di dcf) and dU,' — sin r dr d(j), where 
(j) is azimuthal angle. Obtaining di/dr from SncU's law, 
we finally get 



R' = R 



cosr 



ncosi 

Thus, the modified Fresnel relations are 

RsatEsat,± 2 cosr 



RE_i_ 

RsatEsat,\\ 
REh 



n cos i + cos r 

2 cos ?' 
n cos r + cos i 



The polarization angle is the angle between h and y 
z X X. Thus, 



cosp ^ h ■ y 



1 



smi 



■ (cos i cos adip cos (j)s + sin i sin Ofdip) 



and the sine is found from h x y ~ fcsinp, i.e. 

cos Q!dip sin <ps 



smp ■ 



sin 6 



Under this convention, the angle p is calculated in the 
CCW direction, when viewed from the source of the wave. 

We can choose the polarization angle so that cosp > 0, 
i.e. p € [— 7r/2, 7r/2], by adding to it tt when cosp < 0. 
Then we get 



smp 



cosp 



sin 6 



-sign(cos i cos 9 + sin ctdip) 



1 



sin^ p 



where the argument of sign function has the same sign 
as the previous expression for cosp. 
The electric field components are 

E±^ = Esinp 
E'll = Ecosp 

so that the unit vectors in the directions of E±, and 
k make a right-handed triad. 



3. Polarization and emission angles 

Although E is given by QJ, we need components E± 
and E'li to describe the refraction. Consider a particle 
shower whose direction is described by a dip angle below 
horizon adip and azimuthal angle (in respect to the direc- 
tion toward satellite, calculated clockwise) (f)s (see Fig- 
ure 13). Introduce a coordinate system such that z axis 
is vertical upward, x axis is horizontal in the direction 
of the satellite. Then the unit vector along the shower 
axis is a = {cos adip cos 0s, ~ cos adip sin 0s, ~ sin adip}- 
The unit vector in the direction of emission is fc = 
{sin z, 0, cos?}. The emission angle (between the shower 
axis and the emission direction is found from cos 9 — d- k, 
so that 

cos 9 = cos adip cos 4>s sin i — sin adip cos i 
sin 9 = \/l — cos^ 9 

The polarization angle is the angle between the plane 
containing both d and k and the {x, z) plane. Consider 
h = k X d/\k X d\. Since |A; x a| = sin9, 

h = — — - {cos i cos adip sin chg , 
sm 9 

cos i cos adip cos (j)s + sin i sin adip , 
— sin i cos adip sin 0s } 



4. Antenna response 

The analysis is based on information contained in |64| . 
FORTE satellite has two antennas, A and B, perpendic- 
ular to each other and the nadir direction. Antenna A 
is aligned with the ram (forward) direction. Consider a 
signal coming from azimuthal direction a, and at an an- 
gle b with nadir. Let us choose a coordinate system so 
that z axis is nadir, and the arrival direction is in [x, z) 
plane. Then the antenna directions are given by 

A = {cos a, — sin a, 0} 
B = {sin a, cos a, 0} 

The signal arrival direction constitutes angles a a and as 
with the antennas, which are given by 

cos a A = kA ~ — sin b cos a 
cos ttB = kA — — sin b sin a 

The electric field components parallel to antennas are 

Ea — cos6cosa — sina 
Eb — cos 6 sin a -I- cos a 

We use values Ea/ sina^^i and Eb/ sin as, which are de- 
noted as Ex and Ey in 64], as inputs for the antenna 
radiation diagrams (which are also found in [b^]) to get 
the field recorded by the satellite. 
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